Dust models post-Planck: constraining the far-infrared opacity of dust in the diffuse interstellar medium

We compare the performance of several dust models in reproducing the dust spectral energy distribution (SED) per unit extinction in the diffuse interstellar medium (ISM). We use our results to constrain the variability of the optical properties of big grains in the diffuse ISM, as published by the Planck collaboration. We use two different techniques to compare the predictions of dust models to data from the Planck HFI, IRAS and SDSS surveys. First, we fit the far-infrared emission spectrum to recover the dust extinction and the intensity of the interstellar radiation field (ISRF). Second, we infer the ISRF intensity from the total power emitted by dust per unit extinction, and then predict the emission spectrum. In both cases, we test the ability of the models to reproduce dust emission and extinction at the same time. We identify two issues. Not all models can reproduce the average dust emission per unit extinction: there are differences of up to a factor $\sim2$ between models, and the best accord between model and observation is obtained with the more emissive grains derived from recent laboratory data on silicates and amorphous carbons. All models fail to reproduce the variations in the emission per unit extinction if the only variable parameter is the ISRF intensity: this confirms that the optical properties of dust are indeed variable in the diffuse ISM. Diffuse ISM observations are consistent with a scenario where both ISRF intensity and dust optical properties vary. The ratio of the far-infrared opacity to the $V$ band extinction cross-section presents variations of the order of $\sim20\%$ ($40-50\%$ in extreme cases), while ISRF intensity varies by $\sim30\%$ ($\sim60\%$ in extreme cases). This must be accounted for in future modelling.Comment: A&A, in pres

Physical properties of a very diffuse HI structure at high Galactic latitude

The main goal of this analysis is to present a new method to estimate the physical properties of diffuse cloud of atomic hydrogen observed at high Galactic latitude. This method, based on a comparison of the observations with fractional Brownian motion simulations, uses the statistical properties of the integrated emission, centroid velocity and line width to constrain the physical properties of the 3D density and velocity fields, as well as the average temperature of HI. We applied this method to interpret 21 cm observations obtained with the Green Bank Telescope of a very diffuse HI cloud at high Galactic latitude located in Firback North 1. We first show that the observations cannot be reproduced solely by highly-turbulent CNM type gas and that there is a significant contribution of thermal broadening to the line width observed. To reproduce the profiles one needs to invoke two components with different average temperature and filling factor. We established that, in this very diffuse part of the ISM, 2/3 of the column density is made of WNM and 1/3 of thermally unstable gas (T ~2600 K). The WNM gas is mildly supersonic (~1) and the unstable phase is definitely sub-sonic (~0.3). The density contrast (i.e., the standard deviation relative to the mean of density distribution) of both components is close to 0.8. The filling factor of the WNM is 10 times higher that of the unstable gas, which has a density structure closer to what would be expected for CNM gas. This field contains a signature of CNM type gas at a very low level (N_H ~ 3 x 10^19) which could have been formed by a convergent flow of WNM gas.Comment: 13 pages, 12 figures, accepted for publication in A&

Direct Estimate of Cirrus Noise in Herschel Hi-GAL Images

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.Comment: 4 pages (in journal), 3 figures, Astronomy and Astrophysics, accepted for publication 13 May 201

Statistical properties of dust far-infrared emission

The description of the statistical properties of dust emission gives important constraints on the physics of the interstellar medium but it is also a useful way to estimate the contamination of diffuse interstellar emission in the cases where it is considered a nuisance. The main goals of this analysis of the power spectrum and non-Gaussian properties of 100 micron dust emission are 1) to estimate the power spectrum of interstellar matter density in three dimensions, 2) to review and extend previous estimates of the cirrus noise due to dust emission and 3) to produce simulated dust emission maps that reproduce the observed statistical properties. The main results are the following. 1) The cirrus noise level as a function of brightness has been previously overestimated. It is found to be proportional to instead of ^1.5, where is the local average brightness at 100 micron. This scaling is in accordance with the fact that the brightness fluctuation level observed at a given angular scale on the sky is the sum of fluctuations of increasing amplitude with distance on the line of sight. 2) The spectral index of dust emission at scales between 5 arcmin and 12.5 degrees is =-2.9 on average but shows significant variations over the sky. Bright regions have systematically steeper power spectra than diffuse regions. 3) The skewness and kurtosis of brightness fluctuations is high, indicative of strong non-Gaussianity. 4) Based on our characterization of the 100 micron power spectrum we provide a prescription of the cirrus confusion noise as a function of wavelength and scale. 5) Finally we present a method based on a modification of Gaussian random fields to produce simulations of dust maps which reproduce the power spectrum and non-Gaussian properties of interstellar dust emission.Comment: 13 pages, 13 figures. Accepted for publication in A&

An Imprint of Molecular Cloud Magnetization in the Morphology of the Dust Polarized Emission

We describe a morphological imprint of magnetization found when considering the relative orientation of the magnetic field direction with respect to the density structures in simulated turbulent molecular clouds. This imprint was found using the Histogram of Relative Orientations (HRO): a new technique that utilizes the gradient to characterize the directionality of density and column density structures on multiple scales. We present results of the HRO analysis in three models of molecular clouds in which the initial magnetic field strength is varied, but an identical initial turbulent velocity field is introduced, which subsequently decays. The HRO analysis was applied to the simulated data cubes and mock-observations of the simulations produced by integrating the data cube along particular lines of sight. In the 3D analysis we describe the relative orientation of the magnetic field $\mathbf{B}$ with respect to the density structures, showing that: 1.The magnetic field shows a preferential orientation parallel to most of the density structures in the three simulated cubes. 2.The relative orientation changes from parallel to perpendicular in regions with density over a critical density $n_{T}$ in the highest magnetization case. 3.The change of relative orientation is largest for the highest magnetization and decreases in lower magnetization cases. This change in the relative orientation is also present in the projected maps. In conjunction with simulations HROs can be used to establish a link between the observed morphology in polarization maps and the physics included in simulations of molecular clouds.Comment: (16 pages, 11 figures, submitted to ApJ 05MAR2013, accepted 07JUL2013

ISOCAM observations of the Ursa Major cirrus: Evidence for large abundance variations of small dust grains

We present mid-IR imaging observations of a high Galactic latitude cirrus obtained with the ISO camera ISOCAM at 6" angular resolution. The observations were done with two filters LW2 (5-8.5 microns) and LW3 (12-18 microns) that measure respectively the aromatic carbon bands and the underlying continuum emission from small dust particles. Three 0.05 square degree images sample atomic and molecular sections in the Ursa Major cirrus. These images are compared with Hi, CO and IRAS observations. In such a cloud transparent to stellar light (Av < 0.5) the mid-infrared to 100 microns and the mid-IR emissivity per hydrogen are related to the abundance and the optical properties of small dust particles independently of any modelling of the penetration of the radiation. Within the atomic section of the cloud, the comparison between ISOCAM images and 21 cm interferometric data highlights an enhancement of the mid-IR emitters abundance by a factor ~5 in an Hi filament characterized by a large transverse velocity gradient suggestive of rotation. Furthermore, a drop in the abundance of the same mid-IR emitters is observed at the interface between the atomic and molecular cirrus sections. We propose that these abundance variations of the mid-IR emitters are related to the production of small dust particles by grain shattering in energetic grain-grain collisions generated by turbulent motions within the cirrus and inversely by their disappearance due to coagulation on large grains. At the atomic-molecular interface we also observe a decrease of the Lw2/Lw3 ratio by a factor 2 in the molecular region. This could result from a modification of the dust size distribution or of the intrinsic optical properties of the small dust particles.Comment: 11 pages, 13 figures, better resolution figures to be found in published versio

Magnetic field morphology in nearby molecular clouds as revealed by starlight and submillimetre polarization

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.Comment: 15 pages, 10 figures, submitted to A&

Correlated Anisotropies in the Cosmic Far-Infrared Background Detected by MIPS/Spitzer: Constraint on the Bias

We report the detection of correlated anisotropies in the Cosmic Far-Infrared Background at 160 microns. We measure the power spectrum in the Spitzer/SWIRE Lockman Hole field. It reveals unambiguously a strong excess above cirrus and Poisson contributions, at spatial scales between 5 and 30 arcminutes, interpreted as the signature of infrared galaxy clustering. Using our model of infrared galaxy evolution we derive a linear bias b=1.74 \pm 0.16. It is a factor 2 higher than the bias measured for the local IRAS galaxies. Our model indicates that galaxies dominating the 160 microns correlated anisotropies are at z~1. This implies that infrared galaxies at high redshifts are biased tracers of mass, unlike in the local Universe.Comment: ApJ Letters, in pres